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In Part I of this contribution, a systematic coarse-grained description of the dynamics of a weakly- 
bending semiflexible polymer was developed. Here, we discuss analytical solutions of the established 
deterministic partial integro-differential equation for the spatio-temporal relaxation of the backbone 
tension. For prototypal experimental situations, such as the sudden application or release of a 
strong external pulling force, it is demonstrated that the tensile dynamics reflects the self-affinc 
conformational fluctuation spectrum in a variety of intermediate asymptotic power laws. Detailed 
and explicit analytical predictions for the tension propagation and relaxation and corresponding 
results for common observables, such as the end-to-end distance, are obtained. 

PACS numbers: 87.15.He, 87.15.Aa, 87.16.Ka, 83.10.-y 



I. INTRODUCTION 

Polymer physics has traditionally focused on very flexi- 
ble polymers that admit a highly coarse-grained descrip- 
tion in terms of Gaussian chains and exhibit universal 
physical behavior that can be explained using methods 
from statistical mechanics such as the renormalization 
group and scaling arguments Further research has 
explored new terrain that lies beyond the realm of ap- 
plicability of this highly successful approach either be- 
cause the polymers of interest are too stiff or because they 
are subject to extreme forces. Many of these instances 
have recently appeared in applications involving biopoly- 
mers. Notorious examples are the nonlinear mechanical 
response of DNA Q , which has turned out to be pivotal 
to protein-DNA interactions, and the problem of force 
transduction through the cytoskeleton 0, 01 El j which is 
a major mechanism by which cells explore their environ- 
ment and react to external mechanical stimuli. Clearly, 
in neither of these situations can theorists contend them- 
selves with the convenient Gaussian chain representation, 
but have to resort to more realistic, yet still schematic de- 
scriptions, such as the freely-jointed chain (e.g. for single- 
stranded DNA) or the wormlike chain model (for double- 
stranded DNA, F-actin, microtubules etc.) 

Suspicions that this might entail a substantial loss 
of universality and render systematic analytical ap- 
proaches forbiddingly complex have turned out to be un- 
founded. The wormlikc chain model provides an analyt- 
ically tractable standard model for many of the above 
mentioned new applications, in particular for calculating 
the non-equilibrium dynamical response of stiff and semi- 
flexible but weakly bending polymers to strong external 
fields. As established in Part I of this contribution, the 



weakly bending wormlike chain lends itself to a multiple- 
scale perturbation theory (MSPT) based on a length 
scale separation between longitudinal and transverse dy- 
namic correlation lengths. In the present Part II, we 
demonstrate that the self-affine roughness, acquired by 
the weakly-bending contour in thermal equilibrium, plays 
an analogous role as the more familiar fractal conforma- 
tional correlations in the case of flexible polymers 
The self-similarity of the static conformational fluctua- 
tions entails self-similar dynamics. It manifests itself in 
a variety of intermediate asymptotic dynamic power laws. 
Apart from the restriction to polymers with a (locally) 
rodlikc structure, these predictions are as universal as 
those of classical polymer physics. They are moreover 
derived in a direct way, usually including exact ampli- 
tudes, from a controlled perturbation expansion. 

As a major result of the multiple scale theory devel- 
oped in Part I, we obtained a coarse-grained reformula- 
tion of the free-draining Langevin equations of motion 
of a weakly bending rod in the form of the deterministic 
equation 



d 2 s F(s,t) = -Q (Ag) [F(s,t<t),t] 



(1) 



It describes the long wave-length (all time) dynamics of 
the time-integrated tension 



F(s,t)= / dt'f(s,t') 



(2) 



with £|| being the friction coefficient for longitudinal mo- 
tion and (Ap) [F(s, t < t),t] the average release of con- 
tour length stored in the transverse undulations up to 
time t. Written out in terms of the transverse normal 
mode contributions, the latter reads 



(Atj) (t) 
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where k, t p = n/(k,BT) and C_l are the bending stiff- 
ness, persistence length and friction coefficient for trans- 
verse motion, respectively. The parameter /< = f(t < 
0) =const. allows to take into account a constant pre- 
stress [35| . Upon inserting this dynamical force-extension 
relation into Eq. we arrived at our central result, the 
closed partial integro-differential equation (PIDE) for the 
time-integrated tension F(s,t), 

WM) = (l_ e-^ t+F ^) 

-2q 2 f dte- 2 ' l2 l^ t -^ +F ^- F ^} . (4) 

For convenience, we have made the following choice of 
units: Time and tension, respectively, are rescaled ac- 
cording to 

t -> (±t/K, (5) 
/ - «/ • (6) 

This corresponds to setting k = (± = 1 and ( = 1/2 — . 
As a consequence all variables represent powers of length, 
e.g., t and / are a length 4 and a length -2 , respectively. 

To leading order in the small contour undulations, 
the deterministic coarse-grained tension dynamics, de- 
scribed by Eq. Q, together with the microscopic trans- 
verse equation of motion, represents a valid reformula- 
tion of the constrained Langevin dynamics of a weakly 
bending rod subject to a putative pre-stress. The PIDE 
Eq. J2J is the basis not only for discussing the tension 
dynamics itself, but also the starting point for analyt- 
ical and numerical calculations of the longitudinal and 
transverse nonlinear response of a weakly bending poly- 
mer. It is the purpose of the present Part II to treat the 
former case in detail, while the latter, somewhat more 
complex case is reserved for a future communication Q ■ 
In Sec. m we derive detailed solutions to Eq. for ide- 
alized experimental protocols involving a representative 
selection of external fields. The analytical scaling solu- 
tions obtained for a semi-infinite polymer suddenly pulled 
(or released) at its end, reveal the non-trivial short-time 
phenomenon of tension propagation (Sec. IIII|) . At long 
times, the finite contour length comes in as an additional 
characteristic length scale, which gives rise to additional 
scaling regimes, discussed in Sec. IIVI To make contact 
with experiments, we finally identify the repercussions of 
the tension dynamics on pertinent observables like the 
(projected) end-to-end distance (Sec. 0) and comment 
on novel experimental perspectives brought up by our 
analysis fSec. IVIII . 

II. GENERIC LONGITUDINAL DRIVING 
FORCES 

In general, the tension dynamics depends on how the 
filament is driven externally, i.e., on the boundary and 
initial conditions imposed on Eq. With the defini- 
tion of generic experimental force protocols, this section 



shall provide a framework for the analysis of Eq. 
We introduce the scenarios Pulling, Towing, and Release 
and report on existing investigations. Apart from be- 
ing directly relevant for experiments, we have chosen to 
consider these scenarios because of two properties that 
render them attractive from a theoretical perspective. 
Firstly, they correspond to sudden changes of the en- 
vironment and thus do not introduce an additional delay 
time scale. Secondly, since external forces are assumed 
to act at the ends, these scenarios only change boundary 
conditions and leave the equations of motion unchanged. 
In more complicated scenarios, that involve forces ap- 
plied not only at the ends, these problems show up as 
subproblcms. For instance, if a single point force is ap- 
plied somewhere within the bulk of a polymer, the fila- 
ment can be partitioned into two sections that perceive 
the external force only at their ends. 

1. Pulling 

The polymer is supposed to be free for negative times, 
such that it is equilibrated under zero tension at time 
zero, i.e., we require 

/<=/(«,* <0)=0. (7) 

Then, for positive times the polymer is pulled in 
longitudinal direction at both ends with a constant 
force f j3|| . The corresponding external force density 
f [S(s — L) — S(s)] provides the boundary conditions for 
the tension (see, e.g., the longitudinal equation of motion 
in Part I) 

/(* = 0,t>0) = f , /OM>0) = f . (8) 

Pulling was first considered by Seifert et al. || (SWN). 
They predict that a "large" tension spreads within a time 

t a characteristic length £||(t) = £swN(t) = if/ 2 (ft) 1 / 4, 
from the ends into the bulk of the filament. Their anal- 
ysis neglects bending forces and thermal forces for the 
dynamics, albeit the self-affine thermal initial conditions 
are used (taut-string approximation). The contribution 
of Everaers et al. [3] (EJAM) shed light on the linear 
response to longitudinal forces. Their simulations estab- 
lished a typical propagation length of £\\ (t) = ^ejam(^) = 

£p /i 1 / 8 for weak forces, which was previously predicted 
by Morse H3, and made it plausible by scaling argu- 
ments. Brochard-Wyart et al. (BBG) proposed 
a theory for tension propagation claimed to be valid 
on scales much larger than £ p supposing, however, the 
weakly bending approximation. A quasi-static approxi- 
mation underlies their analysis, in which the polymer is 
at any instant of time considered to be equilibrated with 
the local tension. Applying their results to the situa- 
tion considered here, tension should propagate a distance 

£||(i)«£ B BG(U^ /2 f 3/4 * 1/2 - 

Naturally, the scaling arguments used to predict the 
three different scaling regimes did not address the 
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crossover and the range of validity. Below, we show that 
in fact only two scaling regimes exist. 

2. Towing 

Pulling of a filament can also be studied for time de- 
pendent external forces. A dynamic force protocol of 
particular experimental relevance is given by the con- 
stant velocity ensemble: The polymer is pulled by a time- 
dependent external force f (t) at the left end such that this 
end moves with a constant velocity v [Towing). Besides 
the growth law of the boundary layer, we wish to under- 
stand the time dependence of the external force. We will 
see that both quantities are proportional to each other 
because the external force essentially has to drag a poly- 
mer section of length £\\(t) through the viscous solvent 

with the constant velocity v, hence f (i) ~ £vl\\(t). By 
measuring the time-dependent force in a constant veloc- 
ity experiment one can thus directly monitor £\\(t). Pos- 
sible experimental realizations are outlined in Sec. I VII 

The external force density field corresponding to 
Towing is given by — f (t)S(s) with an external force f (t) 
determined to fulfill the requirement that <9 t r||(0,i) = v. 
Recall from the equations of motion derived in Part I, 
that, up to terms of order 0(e), the gradient of the ten- 
sion is given by the longitudinal friction (as in a rigid 
rod). This implies the boundary condition 

9 s /( S = 0,i>0) = -(v + 0(e) (9) 
f(L,t>0) = 0. 

3. Release 

Release refers to the process "inverse" to Pulling: the 
filament is supposed to be equilibrated at t = under a 
constant pulling force, 

/< = /(M<0) = f >0. (10) 

Then, at t = 0, the external force is suddenly switched off 
and the filament begins to relax. The ends are considered 
to be free for t > 0, 

/( s = 0,t>0) = 0, /OM>0) = 0. (11) 

Release has been discussed by Brochard et al. Ac- 
cording to that work the characteristic size of the bound- 
ary layer, where the tension is appreciably decreased from 
f , should be given by £\\(t) cx ^bbg(^) (the same as for 
Pulling). 

Furthermore, Brochard et al. predict that the tension 
is relaxed as soon as the tension has spread over the whole 
filament yielding a relaxation time t\ for the tension that 

satisfies £ii(i|) = L. This is in conflict with what we will 
find Sec. IIVI where we identify a novel scaling regime of 
homogeneous tension relaxation. 



III. TENSION PROPAGATION 

To unravel the physical implications of Eq. |0J for the 
scenarios introduced above, we begin with the tension 
propagation regime £ii <C L, where the total length L of 
the polymer is irrelevant. In this regime, it is legitimate 
to discuss the dynamics on a (formally) semi-infinite 
arc length interval [0,oo[. Problems like Pulling and 
Release still depend on four independent length scales 
(£ p ,j _1 / 2 ,s,t 1 / 4 ). Yet, it is shown in Sec. UlLAl that. 
Eq. is solved exactly by a tension profile that obeys 
a crossover scaling form depending on only two argu- 
ments, which can be identified as a reduced time and arc 
length variable, respectively. In Scc. lIIIBl we then argue 
that for asymptotically short (<) and long (>) times the 
scaling function reduces to a function of only one scaling 
variable £< = s/£^(t). Our major results concerning ten- 
sion propagation, particularly our classification of tension 
propagation laws £^(t), are summarized in Sec. IIII Cl 

A. Scaling forms 

For each of the generic problems introduced in Scc.lTTl 
we shall see that the tension profile obeys certain 
crossover scaling forms, that cannot be inferred from di- 
mensional analysis. These scaling forms greatly simplify 
the further analysis of the tension dynamics by reducing 
the number of independent parameters. 

To solve the equation of motion, Eq. (@J), for a given 
force protocol, we proceed in the following way. A scaling 
ansatz is postulated and shown to eliminate the param- 
eter dependence in Eq. (@J) and the boundary conditions 
after a suitable choice of length, time and force scales. 
These crossover scales turn out to separate two different 
regimes, a short- and long-time regime, respectively |37j . 

Although being ultimately interested in the tension 
profile f(s,t), it is convenient to first discuss the time 
integrated tension F(s,t), defined in Eq. @, because the 
equation of motion for the tension, Eq. Q), is naturally 
formulated in terms of F(s, t). The physically more intu- 
itive quantity f(s,t) = dtF(s,t) is extracted afterwards 
by a differentiation with respect to time. 

We make the following ansatz for the time integral 
F(s, t) of the tension 

F( M) = f^(A£) , ( 12 ) 

in terms of as yet unknown crossover time and length 
scales if and Sf to be determined below. While a force 
scale f is given explicitely in the case of Pulling and 
Release by the pulling/pre-stretching force, a natural 
force scale 

f = C'VSf (Towing) (13) 

for Towing is provided not unless Sf is fixed. The com- 
bination of variables in Eq. (|13f) represents the force nec- 
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(a, t < 0) = ct = 



(14) 



cssary to drag a polymer section of length Sf (longitudi- 
nally) through the fluid with the imposed towing velocity 
v. 

As long as the polymer is in equilibrium, t < 0, the 
dimensionless scaling function 0(er, r) for the integrated 
tension is zero for both pulling scenarios, but linearly 
increasing with time for Release due to the constant pre- 
stretching force, 

, Pulling / Towing 
t , Release 

The constant c entering the initial condition, Eq. (|14|) . is 
given by c = for Pulling / Towing and c = 1 for Release , 
respectively. 

Since we expect that the signal of a sudden change 
at the end of the polymer, i.e. at a — 0, takes time 
to propagate into the bulk of the polymer, which cor- 
responds to o~ — ► oo, we look for solutions that have a 
time-independent stored length at a — > oo. According to 
Eq. QJ, this corresponds to the boundary condition of a 
vanishing curvature of the tension profile at infinity, 

dl<f>(a -> oo, r > 0) = . (15) 

At the origin, the force, respectively, the gradient of 
the force are prescribed by the considered experimental 
setup, 

r, Pulling 
, Release 
d a (j)(<T = 0, t > 0) = — r , Towing (16b) 

Inserting the scaling ansatz, Eq. Ill2[l . into Eq. Q yields 
after the variable substitutions q — » ?vT > ^ ~~ * T ^f an d 



(ct = 0,t > 0) 



f 3/2/- 1 

dq 
> 



(17) 



3 -2g 2 [g 2 r+0(<T,r)]t f f : 



d ~ e -2g 2 [ 9 2 (t-t )+^(<7,T)-f^(tr,T)]i f f : 



By fixing the scales t$ and st appropriately, 



if = r 



(18a) 
(18b) 



Sf = c -- 1/2 ^/ 2f -i/4 > 
we can eliminate the parameter dependence of Eq. I|17|) 

-,Z7T 



3 



1 



-2<f 



df ( 



-2 9 2 g 2 (T-f) + <^ 



2g 2 [ g 2 r+0(<T,r)]" 

i, (19) 



Note that for Towing, the conditions in Eqs. (|18al I18bj) 
imply the scales 

/ . \-4/5 

^ = \fCh) (20a) 
Sf = (vC 3 e; 2 ) (Towing) , (20b) 



since f depends on Sf via Eq. (|13|l . 



B. Asymptotic scaling 

We have thus removed the parameter dependence of 
the differential equation as well as the boundary and ini- 
tial conditions. The remaining task is to solve Eq. I|19(l 
for 4> under the initial/boundary conditions given by 
Eqs. \ 141 1151 UBI) . and to extract the tension 



f(s,t) 



d t F(s,t) 

s t 

Sf ' tf 
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in terms of the scaling function 

<p(a, t) = d T <j)(a, t) . 



(21a) 
(21b) 

(22) 



Although this remaining task of finding a solution is 
analytically not possible in generality, we now know, at 
least, that it should be a function of only two variables, 
an effective space and time variable a = s/ss and tt , 
respectively. From the conditions in Eqs. i|18al 118 bl) . it 
is seen that the scales Sf and tf are given by a combina- 
tion of the characteristic force scale f (a length" 2 in our 
units) and the persistence length, and hence not simply 
a consequence of dimensional analysis. The significance 
of these nontrivial scales is that they mark a crossover 
in the behavior of the tension. For it turns out that, the 
two-parameter scaling form ip(o~, t) collapses onto one- 
parameter scaling form in the limit of large and small 
arguments, i.e., 



¥>(<t,t) 



a 



for r ^ 1 

> 



(23) 



with a positive and monotonous scaling function </3(£) 
that is bounded as £ — > {0, oo} and exponents a and z 
depending on the problem and the limit, i.e., short or long 
time limit. Equation (|23|l expresses the asymptotic self- 
similarity of the tension profile: by stretching the tension 
profile at a given time in the arclength coordinate a one 
obtains the tension profile at a later time, a property 
inherited from the self-affinc conformational fluctuation 
spectrum of the weakly bending wormlike chain. 

Rewriting the scaling variable in Eq. I|23|) as ct/t 2 = 
s/£»(t) identifies the tension propagation length 



(24) 



In Tab. H] the actual growth laws i\\(t) arc tabulated de- 
pending on the problem and the asymptotic limit. 

Before deriving these growth laws from an asymptotic 
analysis of Eq. H19[). let us give a simple argument as to 
what the exponent z should be on short times. As usual, 
we assume the crossover should occur when the scaling 
variable o~/t z is of order one, hence 



(25) 
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If we further assume, that on very short times the prop- 
agation length £^(t) °f the tension should actually be 
independent of the external force, we can immediately 
infer 



*f(t)oc#V/ 8 



(26) 



which is the correct short time growth law, as will be 
shown in Sec. IIII B II 

The derivations we present in the following are con- 
sistent in the sense that we use assumptions that are a 
posteriori legitimized by the solutions. In particular, we 
exploit Eq. H23|) as a scaling ansatz in order to derive 
asymptotic differential equations for the tension. Those 
equations are then shown to be indeed solved by simi- 
larity solutions of the postulated type. In addition, let 
us assume that the exponent a in Eq. (|23[) is larger than 
-1/2, 



a > -1/2 , 



(27) 



i.e., the tension should increase (decrease) less rapidly 
than t -1 / 2 for r — > (t — » oo). The assumption is 
reasonable for Pulling, because at the ends, we have 
<p(0, r) = 1 and therefore a = in this case. It turns out 
that Eq. (|27fl is correct for all considered problems of ten- 
sion propagation except for sudden temperature changes, 
discussed in Ref. |12| . As a consequence, the approxi- 
mations that are made in the following do not apply to 
sudden changes in persistence length, which is an indica- 
tion that it is an exceptional problem. In Sec. IIVI where 
we consider the scaling regime succeeding tension prop- 
agation, we encounter an asymptotic regime of Release 
characterized by an exponent a = —2/3 as another im- 
portant exception of Eq. i|27|> . 

Since the central PIDE Eq. H19|) is expressed in terms 
of the time integrated tension <fi(cr, r), it is useful for the 
following discussion to reformulate the scaling assump- 
tion, Eq. H23|) . in terms of </>, 



for) 



for r 1 

> 



(28) 



1. Short times (i <C if ) 

In case Eq. (|27fl holds we can linearize Eq. I|19|) for 

short times in 4>. This is at first sight only correct for 
small wave numbers that satisfy q 2 <fi — 0(q 2 T a+1 ) <C 
1. However, upon a closer inspection of the region of 
large wave vectors q > r~( a+1 ^ 2 that do not allow for a 
linearization, it is seen that there is another term in the 
exponent, q 4 T > r" 2 "" 1 > 1 (for t< 1 and a > -1/2), 
which renders the considered exponential essentially zero. 



Therefore, we can approximate Eq. 1191) by 



1 



dq 

J2~^ [q 2 + c 



l-(l-2 9 2 0(a,r))e- 2 ^ 



-2q 2 dt[l- 2q 2 (<t>(a, r) - <f>(a, f ))] e -V(r-f) 
Jo 

' dq ( c 
D 2n \ q 2 (q 2 + c) 



1 - e 



-2q 4 T 



+24>{ct,t) 



1 - 



q 2 + c 



-4<7 / dr <fi(a, f )e 
Jo 



-2g 4 (r-f) 



(29) 



Since r « 1 we can neglect the parameter c in the de- 
nominator of the first term and set the exponential func- 
tion in the second term equal to one. Furthermore, we 
observe that the lower bound of the time-integral can be 
set to — oo by defining 0(er, r < 0) = 0. After the variable 
substitution f — > f + r wc obtain 



%4>(<r,r) 



Q ' <P\P,T) 
q 2 + c 



'dq 



-2<? 4 t _ 1 



2q 



I di(j)(cT,T + r)e 2qif 

J — OO 



(30) 



Now we introduce the Laplace transform of 4> according 
to 



dr ( 



c+ioc 



dz 
2nd 



(a,z) 



(31a) 
(31b) 



such that Eq. (|30|l in Laplace space reads 
d 2 Ma,z) = 



r°°dq 



1 



1 



q 2 (j){a, z) 



q 2 + c 



2q 4 \ z z + 2q 4 
dre~ ZT I df 4>{a,T + T)(2q i )e 2qif 



' dq 



c 



Z(f>(a, z) c(j)(cr, z) 



z(z + 2q 4 ) 2q 4 + z q 2 + c 



= 2~ 3 ' A (z 1 '^ ~ cz- 7 / 4 ) - c 1 ' 2 ^ . (32) 
In the short time limit z>l, this reduces to 

&„<\>(p, z) = 2- 3 / 4 z 1 / 4 . (33) 

Choosing the decaying solution we find for the boundary 
conditions in Eqs. (|15l I16a|) of Pulling 

<f> = <f>p{cr, z) = z- 2 e- a{z/8)1/S (Pulling) . (34) 
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Since Eq. I|32(l is a linear differential equation, we can ex- 
press the solutions for the boundary conditions of Towing 
and Release in terms of (j>p, 

(j) — z~ 2 — 0p(cr, r) {Release) (35a) 

o 

da(f>p(a,T) (Towing). (35b) 



Ultimately, we arc interested in the tension 

ipp(a,T) = <9 T 0_p(er, t) , (36) 

which is given by the inverse Laplace transform of 

z<j) P (a, z), 



<p P (a,T) = 0p[^Js) 

MO = 



ioo+e 7 

az -Hz/8)^ 8 +z 



+f 2iriz 



(37) 
(38) 



After deforming the contour of integration such that it 
encloses the branch cut at the negative real axis, the 
integral in Eq. (|38() becomes 



D dx 8 



sm 



£x sin ■ 



-£x COS ^ 



X 1 - e" 



(39) 



which is to our knowledge not tabulated, but can be eas- 
ily evaluated numerically, see Fig. ^ Upon using the 
known Laplace transform, 



r(i,) 



dTe- ZT T v - x 



9?z/ > , 



(40) 



and Taylor expanding the integrand of Eq. <|38|) one ob- 
tains an expansion of </9_p(£) that is particular useful for 
small £, 



ip P 



(0= E 



(_g-3/8fln 

- n!L(l-n/8) 

/8gN 



(41) 



With an absolute error less than one percent the scaling 
function is approximated by an exponential, 



<£p(0 ~ exp 



23/sr(7/8) 



(42) 



where the pre-factor of £ in the exponent is the initial 
slope d^if\^=o °f the scaling function. 

The tension profiles of the problems under considera- 
tion arc all related to the scaling function </?p(£), 



f(s,t) = ftfp 



hit) 



/(*,*) = f 



1 — <Pp 



(Pulling) (43a) 
(Release) (43b) 



f(s,t)=(v£ ll (t) d£<p P (0 (Towing), (43c) 

Js/l n (t) 



0p(O 




FIG. 1: The scaling function <pp(£) describing the shape of 
the tension profile of Pulling on short times. 



where the boundary layer at time t has the typical size 



£\\{t)=C x 'H x 'H l l s ocl^ MA (t). 



(44) 



The scaling was anticipated in Eq. (|26|) and in an ad 
hoc scaling argument presented in Part I. The result for 
Towing, Eq. I|43c|) . shows that the force at the grafted 
end, at £ = 0, scales like (,v£\\(t). This scaling will be 
shown to hold also in the long-time limit and can be un- 
derstood in the sense that the graft has to balance only 
the drag arising within the boundary layer, since the bulk 
of the filament is not moving longitudinally. Thus, mea- 
suring the force at the grafted end, we can monitor the 
spreading of the tension. This gives special experimen- 
tal relevance to the Towing scenario. In particular, we 
predict 



/(0,t) = C^||W2 3/8 /r(9/8) 



(45) 



for the force at the grafted end on short times, r <C 1. 
The pre-factor in Eq. I|45|) has been found by evaluating 
the remaining integral in Eq. (|43c|) , which yields a Taylor 
series very similar to Eq. (|41|l . 



diMi) = ^ /a E 



(_ 2 -3/8£) 



n+1 



(n+l)!T(l-n/8) 



(46) 



2. Long times (t^> tf ) 

The present subsection deals with the dynamics of the 
tension on a semi-infinite filament on asymptotically long 
times. We identify and interprctc the terms dominat- 
ing the stored length release in this limit. Neglecting 
subdominant terms in the continuity equation, Eq. lfT)l. 
results in differential equations that can be solved by sim- 
ilarity solutions. 

The right hand side of the nondimensionalized PIDE, 
Eq. I|19|) , represents the negative change of stored length 
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in adapted units. It can be written as the sum of two 
terms, A and B, where 



A = 



+ OC 



dq 1 



-oo 2?r + c 



I _ e -2q 2 [q 2 T+<f>(<J,T)} 



(47a) 



B = - f°°^2-2q 2 fdr e -V[9 2 (T-f)+0(<x,T)-<«<x,T)] _ 

(47b) 

We already pointed out in Part I that the term A can 
be interpreted as the "deterministic relaxation" of stored 
length (for the fictitious situation "T=0" , i.e. no thermal 
noise, for t > 0). The term B describes the increase 
in stored length due to the thermal kicks and is strictly 
positive. We analyze both terms separately. 

For pre-stretched initial conditions (c = 1) the long- 
time limit of A follows from setting the exponential to 



.4 



+ 00 



dq 



_ x 2tt q 2 + 1 2 



for c = 1 



(48) 



The term A for r — > oo is nothing but the initially stored 
length, which is completely relaxed after a purely deter- 
ministic relaxation. 

The same reasoning cannot be applied in cases of a 
tension- free initial state (both pulling problems): setting 
the exponential to zero in Eq. i|47a|) for c = yields an 
infrared divergence. The exponential has to be retained 
to render the integrand finite at small wave numbers. 
For the dominant small wave numbers one can, however, 
neglect the term 2g 4 r in the exponent because it is small 
compared to the term 2q 2 $, so that we arrive at the 
asymptotic expression 



.4 



dq 1 
27 ? 



1 



-\Mc, r) 

7T 



3 -2q 2 0( CT ,T) 



for c = . 



(49) 



I.e. the "deterministic relaxation" is dominated by the 
tension term and bending can be neglected (as heuristi- 
cally assumed by Seifert et al. §). A more formal justi- 
fication for Eqs. (|4*51 l4*9)l is given in App. IA II 

The term B, describing the stored length generated by 
the thermal kicks, takes for asymptotically large r> 1 
the form 



B 



T>1 



1 



1 



' dq 

^ q 2 + d T <f>(a, t) ~ 2 v /a T 0(cr, r) 



(50) 



independent of the initial conditions. This is shown 
in App. IA 21 upon using the scaling assumptions in 
Eqs. (13 E3- The result, Eq. gJJ, should not come 
as a surprise. It simply represents the (negative) stored 
length of a stiff polymer equilibrated at the (rescalcd) 
tension d T <f). For a constant force, i.e., d T 4> =const., it is 
obvious that the stored length should saturate for long 



times at the corresponding equilibrium value. But also if 
the tension is varying slowly enough in time (a > —1/2) 
the "noise-generated" stored length can be considered as 
quasi-statically equilibrated with the tension. 

Finally, we combine the "relaxed stored length" ex- 
pressed in A and the "noise-generated stored length" B. 

a. Pulling and Towing: For c = and a > —1/2 
the term A cx \f§ oc r Q / 2 + 1 / 2 is much larger than 
B cx (dr^)- 1 ' 2 = 0(j~ a l 2 ). I.e. the effects of noise 
can be neglected on long times (as presumed by Seifert 
ct al. 8J). In this limit, the thermal noise is merely rel- 
evant in preparing the initial state. The relaxation after 
force application for r 3> 1 is purely mechanical, like for a 
pulled string that is initially prepared with some contour 
roughness (taut-string approximation). If we replace the 
right-hand-side of Eq. i|19fl by the asymptotic form of A, 
Eq. H49|) . we obtain the partial differential equation 



-vw) 



(51) 



for the dynamics of the integrated tension <p. 

Eq. (|51|1 represents a Newtonian equation of motion 
for a particle moving in a conservative force field oc ^fifi 
and can be integrated straightforwardly. In the case 
of Pulling the solution, which satisfies the boundary 
conditions in Eqs. i|15l UGafl of Pulling (in particular 
</>(0,t) = r), is given by 



) = r U/(727tt) 1/4 - 1 



(52) 



for a < (727ft) 1 / 4 and 4> = otherwise. The dimension- 
less tension ip is derived from <f> via differentiation, sec 
Eq. 1221) . and obeys a scaling form, 



f(s,t) = f<p(s,t)=f0 



hit) J ' 



(53) 



with a typical boundary layer size proportional to 

^SWn(^), 

e ll (t) = t 1 ^c 1/2 ^ P / Y /4 (54) 

and a scaling function (p given by 

m= [l-e/(727r) 1 / 4 ] 3 (55) 

for £ < (7271-) 1 / 4 and (p = other wise. Towing starts 
with the same initial conditions (c = 0) but with different 
boundary conditions, Eqs. I|15l I16b|l . Again, we have to 
solve Eq. I|51|l . but now under the boundary condition 
d a 4>(0, t) = — t. The solution is 



T 4/3 

32 



1/3 r 



(187i-r) 



1/3 



- 1 



(56) 



for a < (187TT) 1 / 3 and (f> = otherwise. This implies a 
tension profile of 



f(8,t) = Cvi\\(t)<p 



(57) 
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with the typical boundary layer size 
and the scaling function tp given by 



2tt 



1/3 



I-^/CIStt) 1 ^ 



(58) 



(59) 



for £ < (187T) 1 / 3 and tp = otherwise. As on short times, 
the absolute value of the reduced tension at the left end is 
proportional to the size of the boundary layer. Its precise 
value is predicted to be 



/(0,i) = (2 7 r/3) 1/3 C^||W(xi 1 / 3 , 



(60) 



and should be directly accessible to single molecule ex- 
periments. 

b. Release: The stored length release is asymptoti- 
cally given by 



(Ag) ocA + B 



1 



1 



(61) 



This expression for the change in stored length can also 
be directly obtained if one assumes that the filament was 
at any time equilibrated with the current tension tp = d T <p 
{quasi- static approximation). Our derivations show, that 
this assumption, used by Brochard et al. [Tl| . is only valid 
in the long time limit for the scenario of Release. The 
dynamics of the integrated tension is then described by 



01 



2 2y/d r <K*,T) ' 

or, in terms of the tension tp = d T <j){a, r), 
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(62) 



(63) 



The solution satisfying the correct boundary conditions, 
Eqs. (|15l I16a(l , is given by a scaling form 



<p(s,t) = tp 



(64) 



with the typical boundary layer size now growing like 

^|,(t)=t 1 / a r 1 / 2 ^ a f 3 /*oci B BG(t) (65) 

and a scaling function tp(£) satisfying the ordinary differ- 
ential equation 



dl0 = -l&- 3/ %<P 



(66) 



The scaling function depicted in Fig. [21 was already ob- 
tained numerically in Ref. The slope at the origin 
is c\<^ =0 w 0.6193. 



0.1 



FIG. 2: The scaling function <p(£) for Release obtained by a 
numerical solution of Eq. (I66H . 



C. Tension propagation (summary) 

This section summarizes the picture of tension prop- 
agation that emerges from the above solutions of the 
dynamical equation for the tension, Eq. for sudden 
changes in boundary conditions. 

For the considered problems Pulling, Towing and 
Release, we have shown that the tension profile in units 
of f has to obey a crossover scaling form tp(<j, r) depend- 
ing on a reduced time variable r = t/tj and a reduced 
arc length variable a = s/sr. The scaling function tp 
describes how sudden changes of the tension at the ends 
spread into the bulk of the polymer. Its crossover struc- 
ture and the expressions for tj and Sf ~ £\\{tf ) are consis- 
tent with our heuristic discussion of Pulling in Part I. In 
the limits t< 1 and r ^> 1 the function tp(o~, r) assumes 
a simple (one-variable) scaling form 
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Mr* 



(67) 



The notation ^ indicates that the asymptotic form of tp, 
a, and z will generally not only depend on the kind of ex- 
ternal perturbation applied, but will also differ for times 
t ^ tf . Rewriting <j/t z = s/i\\ identifies the tension 
propagation length £ti=Sf t z . 

For t <C tf Eq. Q could be linearized in / and the 
scaling function (^ < was obtained analytically. Whereas 
<p < depends on the considered force protocol, the cor- 
responding exponent z < = 1/8 is independent of the 
boundary conditions. As established by our heuristic dis- 
cussion of Pulling, this is due to the relaxation of modes 
with Euler forces l^ 2 f much larger than the external 
force, for which the equilibrium mode spectrum is hardly 
perturbed by the external force. The self-affinity of the 
equilibrium mode spectrum translates into a self-similar 
relaxation dynamics. The dynamic exponent z for the 
growth of the boundary layer could already be antici- 
pated from requiring < y 5 < to become f —independent as 
in linear response, see Eq. (|26|l . The short-time dynam- 
ics for strong external force is thus closely related to the 
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TABLE I: Asymptotic growth laws for the dynamic size i\\ (t) 
of the tension boundary layer. 



Problem 


t < t f 


tf <t<4 


Pulling 






Towing 






Release 




Vvcf 8/4 * va 



linear response. Note, however, that the limit f — > is 
problematic, as it does not interchange with e — ► 0. Our 
identification of arc length averages with (local) ensemble 
averages in Part I breaks down for f < (£/ '£ p ) 1 'H~'' 'Z 16 , 
where fluctuations in the tension become comparable to 
its average value. In fact, extending Eq. Q to linear re- 
sponse amounts to an uncontrolled factorization approx- 
imation (fr 2 ) — * (/) ( r i)- Even in the stiff limit the 
linear longitudinal dynamic response remains an open 
problem. The limit where fluctuations in the tension be- 
come important and its consequences will be detailed in 
Sec. EH 

For t 3> U the dynamics becomes nonlinear in the 
external force and starts to depend on the kind of ex- 
ternal perturbation and on how precisely it is applied 
to the polymer. Previously predicted power laws were 
recovered from Eq. Q by employing different approx- 
imations to its right hand side. In the taut-string ap- 
proximation of Ref. Q, one neglects for t > bending 
and thermal forces against the tension, i.e., one drops the 
q 4 — term in the relaxation time r q = q 4 + fq 2 of a mode 
with wave number q and sets £ p — > oo for positive times 
(i.e. — 0). The complementary quasi-static approxi- 
mation of Ref. amounts to the omission of memory 
effects, i.e., to the assumption of instantaneous equili- 
bration of tension and stored length (as it would be the 
case for vanishing transverse friction, £j_ — > 0). In cases, 
where either of these approximations applies, a power- 
law dispersion relation combines with a sclf-affinc mode 
spectrum to produce self-similar tension dynamics. Our 
analysis of Eq. Q showed that either of these approxi- 
mations becomes rigorous in the intermediate asymptotic 
regime defined by t 3> tf , £.\\ <C L. The quasi-equilibrium 
approximation applies to Release and the taut-string ap- 
proximation to Pulling. We could rule out the applicabil- 
ity of the taut-string approximation for Release and of the 
quasi-static approximation for Pulling [Tlj and Towing. 
The "pure" scenarios of self-similar dynamics are sum- 
marized in Tab. ^ and Fig. EI 



IV. TERMINAL STRESS RELAXATION 

Up to now, we have considered the growth of the 
boundary layer in a stiff polymer that has a (formally) 
semi- infinite arc length parameter space, s € [0,oo[, 




FIG. 3: Schematic of the tension propagation laws £« (t) oc t z 
on a double-logarithmic scale. At tf = f ~ 2 they crossover 
from a universal short-time regime to (problem-specific) 
tension-dominated intermediate asymptotics, except for weak 
forces, f < f v /L*. The propagation ends when £y (t) w L. 



which is an idealization. However, the foregoing discus- 
sion obviously applies to a polymer of finite length L for 
sufficiently short times: As long as the size of the bound- 
ary layer is much smaller than the total length L the 
presence of a second end is irrelevant to the boundary 
layer at the first end. The time where the boundary lay- 
ers span the whole polymer marks the crossover to a new 
behavior. For definiteness, we define the crossover time 

4 by 

e\\(t[) = L. (68) 

What happens for t > t\l The straightforward way 
to answer this question is to solve for the intermediate 
asymptotics of the PIDE, Eq. (0J, for a polymer of fi- 
nite length. The finitcness of L amounts to replacing 
the boundary condition d 2 f(s — > oo) = by the correct 
problem-specific one, i.e., by 

f(L,t>0) = f Pulling 
f(L,t>0) — Release, Towing. 

One could now proceed as in Sec. IIII Bl by identify- 
ing proper scaling forms and extracting their asymp- 
totic behavior. As compared to the semi-infinite poly- 
mer limit, this procedure is more complicated for a fi- 
nite polymer because of the additional scaling variable 
£±(t)/L. Therefore, we prefer to take the following short- 
cut, which consists in two steps. 



1. Trivial tension profiles for t > t,. 

In the heuristic analysis of Part I, we found that or- 
dinary perturbation theory (OPT) should become valid 
for times larger than some crossover time i*. With the 
exception of Towing, which we discuss separately below, 
the introduced scenarios treat both ends equally, such 
that the tension profile in OPT has a trivial time and 
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arc length dependence, namely 

/° PT (M) =/(0,t) = const. , (69) 

up to small terms of order 0(e). For constant tension we 
can easily extract the longitudinal dynamics from Eq. 
The corresponding predictions for the evolution of the 
end-to-end distance will be discussed in Sec. IV Bl 



2. Possibility of homogeneous tension relaxation for 
t\ <t < U. 

Up to now, we have argued that the tension propa- 
gates for t <C t\ and is constant in time and space for 
( > ( t . There remains the question whether there is a 
non-trivial regime of homogeneous tension relaxation in 
the time-interval jtt; i*]. Using the systematic approach 
outlined in Add. lO to determine t*, we actually find for 
most of the problems that t* — V L , i.e., there is no scaling 
regime between tension propagation and the stationary 
tension profiles dictated by OPT. The i?e/ease-scenario, 
however, provides an important exception, as it allows 

for a time scale separation v L <C t*, as we demonstrate 
cxplicitely below. For intermediate times the tension re- 
laxation is shown to exhibit a novel behavior with an 
almost parabolic tension profile and an amplitude that 
slowly decays in time according to a power law. 



A. Release for large pre-stretching force 

Let us first determine the time at which OPT be- 
comes valid, for Release. In the OPT regime, the tension 
should be so small, that we can calculate the change in 
stored length accurately by means of Eq. © (with the 
pre-stretched initial conditions of Release, i.e. /< = f) 
under the assumption of a vanishing tension, / OPT = 0, 




Here, we have substituted q — > qt l / A and replaced f by 
U ~ 1//2 in order to obtain Eq. JTOdJ). In the final step, 



we took the long-time limit t 3> U ■ Let us now de- 
termine the order of magnitude of the variation <5/ of 
the tension due to the longitudinal friction. Estimating 
Sf ~ C (Ag) /(tL 2 ) from Eq. (JTJ), we obtain 

5f*&-^. (71) 

As in the heuristic discussion of Part I, we find a diverging 
tension in the limit t — > for the OPT result, so that the 
OPT result can only be valid after the effect of Sf on the 
evolution of (Ag) (t) can be neglected. For the particular 
case of Release, this time can be determined as follows. 
As discussed in Part I, the stress-free dynamics is at the 
time t characterized by relaxation of modes with wave 
length £±(t) ~ i 1/4 . When 5f is larger than the critical 
Euler Buckling force £±(t)~ 2 corresponding to the length 
£±{t) we expect that the tension cannot be neglected for 
the evolution of (Ag). Hence, the above result (A~g) oc 
Lt^ 4 /e p , obtained fr om OPT, can only be valid if 

^(t) 2 <5/c^t- 1 /4 <<l! (?2) 

i.e., for long enough times 

t » U = L 8 /£ 4 p . (73) 
The time t*(L) is obviously not identical with the time 

tl = L 2 £~ 1 ir m (*»*f) (74) 

it takes for the boundary layer to spread over the fil- 
ament. To compare them, we first notice that tl ^> 
tj = f ~ 2 implies that the polymer must have been pre- 
stretched by a large enough force, 

f > £p/L 4 = e- 2 L- 2 , (75) 

larger than the pre-stretching force at least necessary to 
enter the propagation regime tj <C t <C! if,. By using the 
estimate in Eq. (|75() we can compare i* and 

= f (f) 3 »fr 3/2 -i- (76) 

It is seen that the time window i| <C t <C i+ grows with 
the pre-stretching force, which means, in particular, that 
it describes the limit of an initially straight polymer. 

To determine the physics of this novel regime, we have 
to solve the equation of motion for the tension on the 
finite arc length interval [0; L] using the approximations 
developed in Sec. IIIIBI There, we found that in the 
limit t 3> tf the tension profile of Release is described by 
Eq. JBSJl 

dl V =\^ 2 d rV . (77) 
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The right-hand-side represents the time-derivative of the 
stored length in the quasi-static approximation. Going 
back to variables /, s and t, Eq. (|77jl takes the form 

o 2 j = \j-r 3/2 d t f. (78) 

For the following, we assume that this quasi-static ap- 
proximation is not only valid in the tension propagation 
regime for tj <C t <C i|, but also for longer times until 
the OPT regime begins, if C t « t t . This is justified a 
posteriori in App. [Bj where it is shown that the change 
in stored length for the solution /(s, t) of Eq. I|78|) can in- 
deed be calculated quasi-statically for times t <C t+. We 
solve Eq. (|78|l for the boundary conditions 

f(s = 0,t) = 0, f(s = L,t) = (79) 

and the initial conditions 

/(s,0) = f , for < s < L . (80) 

In Sec. 1111 li"2l we solved the differential equation Eq. I)77|) 
with the scaling ansatz / = f tp[s/£\\ (t)] numerically for 
the correct initial condition, but ignored the boundary 
conditions of at one end by sending L — ► oo. In contrast, 
we now look for a simple solution that obeys the tension 
boundary conditions exactly, at the expense of a possible 
mismatch with the initial conditions. To this end, we 
make the product ansatz 

f(s,t)=g(t)h(s) . (81) 

Separation of variables yields 

^g- 5/2 dtg = c = Vhh" . (82) 

Choosing C = — 1/(6L 2 ), for convenience, we find the 
long time asymptotics 

m ~ {0) ' • (83) 

The spatial part obeys a scaling form h(s) = h(s/L), 
where h satisfies the equation 

dfh(0 = ~\h- 1/2 (84) 
with the boundary condition 

h(0) = h(l) = . (85) 

The analytical solution can be found in a standard way. 
The main characteristics are the slope at £ = 0, 

%/i(£)|f = o = 12~ 1/3 w 0.4368 . (86) 




FIG. 4: The initial (t <C t*) time evolution of Release. A 
regime of (slow) tension relaxation begins after the sudden 
change in boundary condition has propagated through the 
filament. 

and the maximum value of h, 

h(l/2) = ± (J\ ' w 0.0819. (87) 

Up to now we have investigated the quasi-static ap- 
proximation in the tension propagation regime [t <C t\) 

and in the regime of tension relaxation (t 3> t L ) sepa- 
rately. In order to illustrate the crossover, we have also 
solved the corresponding PIDE, Eq. H78(l . numerically. 
The result shown in Fig. IIV Al unveils the transient na- 
ture of the tension propagation regime. 



V. NONLINEAR RESPONSE OF THE 
PROJECTED LENGTH 

After having discussed the rich tension dynamics of 
stiff polymers in detail, we wish to derive its consequences 
for pertinent observables in order to make contact with 
experiments. Although the tension may in some situ- 
ations be monitored directly, see Sec. IVII a more con- 
ventional observable is the longitudinal extension R\\ (t) 
of the polymer, which is defined to be the end-to-end 
distance projected onto the longitudinal axis. Tension 
dynamics strongly affects the nonlinear response of the 
projected length, which shall be detailed in the following 
for the force protocols Pulling, Towing and Release. 

The average temporal change in the projected length 
R\\ (t) is directly related to the stored length release, 

<Ai?||) (*) = -/ ds (Ag) (a, t) + o(e) , (88) 
Jo 
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which was already noted in Part I. As long as modes with 
wave length on the order of the total length are irrelevant 
(i.e. both ends of the polymer are not correlated), 



t±(t) « L , 



(89) 



half-space solutions for (Ag) may be used to evaluate 
Eq. JSH||. 

Recall from Part I, that (Ag) can be decomposed into 
a bulk contribution (A7j) and a term influenced by the 
boundary conditions, which vanishes under a spatial av- 
erage. We show in Sec. IV II that in the cases of hinged 
and clamped ends, as opposed to free ends, the contri- 
bution of the boundary term to the integral in Eq. I|88l) 
is subdominant in the limit t —> 0. Nevertheless, these 
boundary effects represent important corrections that 
should be taken into account in any experimental sit- 
uation (i.e. with finite t) of Pulling and Towing. 

At first, however, let us consider the bulk contribution 
to the end-to-end distance, 

(AS,,) (t) = [ ds (Ag) [F(s, t<t),t] , (90) 
Jo 

which is universal in the sense that it does not depend 
on the boundary conditions of r_i_. In Eq. I|9U|) . and in 
the following, we neglect contributions of higher order 
than Ag = 0(e) to the projected end-to-end distance. 
According to Eqs. I|9UI IT|). the bulk change of the end- 
to-end distance is given by the arclength integral of the 
curvature of the time integrated tension F(s,t), 



(AR ll )(t)=C 1 / dsF"(s,t) 



= C l [F'(L,t)-F'(0,t)} 



(91a) 
(91b) 



The predictions are presented separately for the regime 
of tension propagation (t <C if,) and the regime of OPT 
(t where the tension profiles are flat and time- 

independent up to sub-leading terms. For weak forces, an 
explicit expression for times t < t* is given in Sec. IV Dl 
which captures the crossover from tension propagation 
to the tension saturated OPT regime. As discussed 
in Sec. IIV Al Release turns out to also have an addi- 
tional time- window t\ < t < t*, whose consequences for 
the end-to-end distance is described in Sec. IV CI Since 
Towing is the only problem, in which both ends do not 
behave in the same way, we discuss this "asymmetrical" 
problem separately in Sec. IV Gl 



A. Tension propagation regime (( < f[) 

In the time domain of tension propagation, where the 
boundary layers are growing in from both ends (t <C t\), 
we may use the tension profiles for the semi-infinite 
(pseudo-) polymer from Sec. IIII Al (labeled by F^ here) 



TABLE II: Universal bulk contribution (AR\\) (t) to the dy- 
namic change of the end-to-end distance in the limit t <C t\ <C 



Problem 



Pulling 



Release 



t < tf 




T(15/8) 



t > t f 



1/4 f 3/4 



V 71 - 




-2.477 



to approximate Eq. (|91b|) by 



<AE,|)(t) *~ L -2r^(0,t) 



(92) 



for scenarios with two equally treated ends. Now, if t 
falls into a regime where the tension exhibits scaling, 



(93) 



we immediately obtain from Eq. I|92|l the power-law 

(AE|, ) (f ) cx -t a+1 - z 2d£ Foo (£ = 0) (94) 

for the growth of the end-to-end distance. The pre- 
factors can be calculated for all cases, because the scal- 
ing functions are known. In this way we obtain the list 
Tab. ITU of growth laws. 

B. OPT regime (t > t*) 

When calculating the release of stored length in OPT 
the tension profile / OPT = f(s = 0) = const, is assumed 
to be stationary and flat, so that 



(Ai?n)W = -L(Ae)[/ OPT ,t] , 
for t > t* . 



(95) 



In the case of Release, the quantity (Ag) (f OPT = 0, t) 
has been explicitely calculated in Eq. I|70() . For Pulling, 
(Ag) (f OPT = f,t) can be evaluated from Eq. © in a 
similar straightforward manner, because the tension is 
spatially constant. The corresponding growth laws are 
summarized in Tab. Hill 

Compared to the tension propagation regime, the 
growth laws in the OPT regime are slowed down, sec 
Tabs. |nj IIHI Actually, for all cases except Release the 
corresponding growth laws obey 

/ a -f: \ OPT / . — > MSPT 

(AE,|) oc(Ai?,|) r* . 

This can be understood in terms of the scaling arguments 
used Part I. There, we took tension propagation hcuris- 
tically into account by assuming that the stored length 
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TABLE III: Universal bulk contribution (AR\\) (t) to the dy- 
namic change of the end-to-end distance in the limit t* <C 

t<ztt. 



Problem 



Pulling 
Release 



t < t 



if 



if (t 
£ P F(7/4) \2 



3/4 



t > t 



2ft 



r(i/4) £ P 



release, as given by OPT, is restricted to the boundary 
layer of size £\\(t). Then, one has 



OPT 



<Ai? it p - ~ -*„(<) (Ae) (/ Ui " r ,t) for * « ^ (96) 
as compared to 

(Ai?n) MSPT ^-L(Ae)(/ OPT ,t) fart»t*. (97) 
This conforms with the heuristic rule noticed above, 



,a+l-2z 



(98) 



C. Release in the limit f" C t < t t 

The intuitive rule in Eq. I|98|) fails for Release in the 
limit t 3> if , which indicates that this is an exceptional 
scenario. As discussed in Sec. IIV Al = L 8 /£* cannot 

be identified with t{ = Z, 2 £" 1 Cf" 3/2 in this case. There 

exists a time window t\ <C i <C f* that expands in the 
limit of large forces, f ^> £ p /L 4 . We have shown that 
the tension exhibits homogeneous relaxation in this novel 
regime. With the slope of the tension at the ends 
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the growth law 

(ASiiXtJ^-lS 1 ^ 



Lt 

M, 



2/3 



l I e p t 



2/3 



(99) 



1/3 



(Release) (100) 



follows from Eq. ij91b(l . We expect the growth law 
(Ai?||) cx t 1 / 3 during homogeneous tension relaxation to 
hold even for chains with L t p . The example of re- 
tracting DNA will be discussed as an experimental out- 
look in Sec. IVII The exponent 1/3 coincides with that 
obtained by an adiabatic application of the stationary 
force-extension relation to a "frictionless" [l^ poly- 
mer with attached beads at its ends |l7j . 

For times t > f* = L 8 /£ 4 the growth law in Eq. (fTUUf) 
crosses over to the one noted in Tab. IIIII Interestingly, 



both growth laws appearing for i > f[ are independent of 
the initial tension f . In both cases, the initial conditions 
are completely "forgotten" once the tension has propa- 
gated through the whole polymer. An overview over the 
time scales separating the diverse regimes for Release (as 
compared to Pulling) will be given in Scc. lVFl 



D. Pulling and Release for small forces 

Provided the external force is smaller than the criti- 
cal Euler buckling force of the polymer, f < L~ 2 , the 
crossover time U exceeds the terminal relaxation time 
if, hence the linearized PIDE, Eq. (|32[) . applies through- 
out the contour relaxation. The linearity allows to solve 
the PIDE for a polymer of finite length, and we ob- 
tain an analytic description of the crossover between the 
asymptotic power laws AiJii (t) oc i 7 / 8 for t -C t" L and 
AR\\(t) oc i 3 / 4 for t* <i < ti (cf. Tabs. El EH. To this 
end, let us first reinstall original units (which are bet- 
ter adapted for the present purpose) into the linearized 
PIDE, Eq. (|22J>, 



dlF{ Sl z) 



2 3 / 4 4, 



F(s,z) , 



(101) 



where F(s, z) is the Laplace transformation of the inte- 
grated tension, 



F(s,z) 



dte~ zt F{s,t) . 



(102) 



For the boundary conditions of Pulling, Eqs. (|15l I16afl . 
the solution to Eq. I|101|) is given by 



F(s,z) = fz- 



cosh 


'pt p )W{z/8)V*(s-L/2j 


cosh 


\t/£ p y/2( z /&y/8L/2_ 





(103) 



This implies a growth of the end-to-end distance of 
(AR ll )(z) = -2C 1 d s F\ s = 

2 5/8 f 



Z lf V8((£ p )l/2 



■ tanh 



1/2 



z\V s L 
8/ 2 



(104) 



The Laplace back-transform of Eq. (|105|l takes the form 
r-ioo+e dz 



<AJJ,|)(f) 



fC 3 L 7 



2m 
T 



F e-*<AH||>(t) 



(105) 



where T(r) is a scaling function, given by 

XT ^ 
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T(r) = - — / ds 
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(106) 
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FIG. 5: The scaling function T(r) in a double logarithmic 
plot. The crossover of the scaling behavior is close to the 
time t c = 6.55 x 10 -4 where the asymptotic short-(blue) and 
longtime (red) asymptotics cross. 

We have depicted T(r) in a double logarithmic plot in 
Pig. (Sj together with the corresponding asymptotics from 
Tabs. HTl Iffll Note that the time t c sa 6.55 x 1CT 4 where 
the asymptotic lines cross is a good indication for the 
crossover occurring in the exact solution. Assuming that 
this also holds quite generally, it is possible to obtain esti- 
mations of when the crossover between linear and nonlin- 
ear regimes should occur in the natural time units tf . To 
this end, one simple equates the corresponding asymp- 
totic power laws (including the exact pre-factors) for the 
growth of the end-to-end distance. As for the present case 
with r c sw 6.55 x 10 , these crossover times are typically 
not of order one in natural units because of the numer- 
ical proximity of the exponents of the asymptotic power 
laws, which are 7/8 and 3/4 in the present case. In a 
given experimental situation, one should therefore check 
carefully which regime is expected by comparing exper- 
imental time scales with these unusual crossover times. 



Tension fluctuations exceed the applied force in magni- 
tude at any given time for small enough external forces. 
Hence, in the limit f — > while e <g; 1 is fixed, the tension 
cannot be considered as a deterministic quantity. Recall 
however, that our MSPT analysis in Part I was based on 
the limit e — > while f is fixed, and only in this limit the 
self-averaging argument of Part I applies. Extending our 
results to the usual linear response limit corresponds to 
the uncontrolled approximation fr'± — ► f (r'±). To ana- 
lyze this limit more carefully, one has to solve the stochas- 
tic PIDE that obtained in Part I before the self-averaging 
argument was employed. However, since the 7/8— scaling 
of the fluctuations has already been confirmed in simu- 
lations 0, we expect that such a more rigorous analysis 
would yield the same scaling but a pre-factor different 
from the one of the deterministic short-time law (Tab. lTT|l . 

F. Release versus Pulling (overview) 

The diverse regimes and their range of validity are sum- 
marized in Fig. El for the Pulling and the complementary 
Release problem. It depicts the crossover time scales as 
a function of the externally applied tension f . The line 
U ~ f ~ 2 separates "linear" from "nonlinear" behavior. 
The symmetry of the graph for t < tt with respect to the 
f = 0-axis indicates that the scenario of Pulling can in- 
deed be considered as the inverse scenario of Release for 
weak forces / <C f c , or, more generally, on short times. 
This symmetry is lost in the nonlinear regime. The grow- 
ing importance of uniform tension relaxation for Release 
with increasing initial tension f becomes particularly ap- 
parent on the log-scale of the figure. Due to the particular 
choice of units (t/L 4 and f / L~ 2 ) the regimes of nonlinear 
growth of the boundary layers appear relatively narrow in 
Fig. El Which of the various regimes will pre-dominantly 
be observed in measurements actually depends strongly 
on the ratio L/£ p and on the experimentally accessible 
time scales. 



E. FDT and linear response 

According to the fluctuation-dissipation theorem 
(FDT [3 111), the fluctuations in the end-to-end dis- 
tance should be related to the linear response of the 
polymer. It is tempting to interpret the linear short- 
time regime, where (Ai?) ~ f i 7/,8 (C^ p ) -1 / 2 , as linear re- 
sponse in the usual sense, in which case the fluctuations 
of the end-to-end distance should, according to the FDT, 
scale as (^AR 2 ^ ~ t 7 / 8 £ -1 / 2 ^p 3 ^ 2 in equilibrium. How- 
ever, fluctuations of this strength are inconsistent with 
a deterministic tension dynamics, because they generate 
friction forces over a scale £||(t) ~ t 7 / 8 (VC) 1/2 and thus 
imply tension fluctuations of magnitude 

Sf ~ Ce {] (i) (AR) /t ~ (C/£ P ) 1/4 t- 7/w (107) 



G. Towing 

Towing was excluded from the preceeding discussion 
because it is the only "asymmetric" problem, as both 
ends of the polymer behave differently. In the tension 
propagation regime (i <C i|), the left end is constrained 
to move with constant velocity, while the right end expe- 
riences no driving force. Hence, we have 

(Aii||)(i) =vt , fort<4- ( 108 ) 

After the boundary layer of non-zero tension has reached 
the free end, t 3> i|, the right end starts to move. Then, 
the tension profile becomes linear as for a straight rod 
dragged through a viscous solvent. The further contour 
relaxation is up to pre-factors identical to the Pulling 
problem for t t| w t*. 
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FIG. 6: Characteristic times (logarithmic scale) for Pulling 
and Release against the applied external force f (linear scale). 
The time £* (stars) separates regions where ordinary pertur- 
bation theory (OPT) applies (dark shaded) from regions (light 
shaded) of linear (hatched) and nonlinear tension propagation 
and from homogeneous tension relaxation (white). Whereas 
longitudinal friction is negligible for t > i*, it limits the dy- 
namics for t < t*. The innermost funnel-shaped region indi- 
cates the regime where the tension fluctuations are important 
(defined by / < Sf with Sf given by Eq. fiE) ). 



H. Pushing 

For completeness, we mention the scenario of a fila- 
ment being compressed by external longitudinal forces. 
This scenario has some subtleties. Pushing increases 
the stored length exponentially for t 3> if by virtue of 
the Euler Buckling instability and generates a situation 
where the weakly bending approximation is not valid any- 
more. Then, hairpins are generated |2dj| and for t tf 
the rigidly oriented driving forces pull on those hair pins. 
Our theory is only applicable at short times. For t<t/ 
the response of the system is linear in the driving force, 
irrespective of the sign. 



I. Boundary effects 

Up to now, we only discussed the change (Ai?||) (t) 
in projected end-to-end distance corresponding to the 
change (A'g) of the stored length in the bulk. For 
a hinged, respectively, clamped semi-infinite polymer, 
the as yet missing boundary contribution X h ' c (t) = 

(Ai?J /c )-(Ai?||) is given by 



X h ' c (t) = ±2 ds 

J0 JQ " L P 



oo />oo 



dq \ e -2<? 2 [<? 2 t+^(M)] _ i 



q 2 + /< 



+ 2q 2 [ dte- 2q2 ^^ +F ^- F ^^\co S (2q S ) . 



TABLE IV: Boundary contribution to the change 
(^AR^^ (i) of the end-to-end distance for hinged, re- 
spectively, clamped ends in the limit tj; <C t. 



Problem 


(AR 1 '/ ) (t) - (AT?,,) (t) 


Pulling 


±2tf/4 


Towing 


±2tf/£ p 


Release 


( ±2tf /t p , for /<«L" 2 


\ 0, for /<»L- 2 



As discussed in Part I, the boundary dependent term of 
the stored length decays on a length scale of 0(1) due 
to the cosine factor. Since the tension decays on a much 
larger length scale of order ©(e -1 / 2 ), it is permissible 
to use the (integrated) tension at s ~ to evaluate the 
arclcngth integral in Eq. 1(109(1 . 

Upon replacing F(s,t) — > F(0,t) and using 



ds cos(2qs) = (n/2)S(q) , 
the integral in Eq. 1)109(1 can be evaluated, 

g^O £ p q 2 + /< 
which vanishes, unless 



./'< 







x h / c {t) = ± 



2F(0,t) 



(110) 



fill) 



(109) 



On the semi-infinite arclength interval, these boundary 
contributions are therefore nonzero only for Pulling and 
Towing, and are summarized in Tab. II Vl Note that, the 
boundary effects of clamped (hinged) ends tend to reduce 
(increase) the longitudinal response of the polymer in 
comparison to the bulk response. This may be explained 
as follows. Close to a clamped end, a polymer is more 
stretched out than in the bulk because r'f_ is constraint 
to approach zero at the end. As a consequence, the end 
portion of the polymer is less able to store or release 
excess length. For hinged ends, the boundary conditions 
act just in the reverse direction. 

The strict vanishing of the boundary term for any fi- 
nite /< and the discontinuity at /< = is a consequence 
of the assumed infinite half-space. For a polymer of finite 
length, it turns out that the boundary term approaches 
zero for pre-stretching forces larger than the critical 
Euler-Buckling force, /< ^> f c = L~ 2 . This can be 
seen by studying the finite integral J Q L ds ((Ag h / C ) (s, t) — 
(Af>) (s, t)). For forces /< / c , the integrand satu- 
rates at a plateau of magnitude 0(f^/ 2 t/£ p ) for £±(t) < 

— 1/2 

s < f < before it finally decays to zero. Within the 
semi-infinite integral, the integral over this long plateau 
cancels the contribution stemming from s < £±(t), as 

required by Eq. (|110fl . However, for ^ 2 the con- 

tribution from the plateau may be neglected. As a con- 
sequence, the value of the integral is for /< <C f c given 
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by ±2f < t/£ p . This asymptotic behavior is important for 
Release, as noted in Tab. UVl 

Upon comparing Tab. II VI with Tabs. ITT1 IIIII one may 

think that boundary effects are always subdominant in 
the short time limit. However, our calculations were spe- 
cialized to hinged or clamped boundary conditions. In 
many experimental situations, one has to deal with free 
boundary conditions. Somewhat tedious but rigorously, 
these boundary conditions can be taken into account by 
means of the correct susceptibility, which can only be 
given in terms of an integral. Here, we discuss effects 
related to free boundary conditions on a heuristic basis 
and show that they generate a dominant contribution to 
the change in the end-to-end distance for Pulling, which 
is proportional to £ 3 / 4 . 

The argument is based on the observation, that exter- 
nal forces that act in the longitudinal direction (||) while 
the polymer is free, automatically introduce small (of or- 
der (Die 1 / 2 )) transverse forces at the ends. This follows 
from the force balance equation of a semiflexible polymer 
(cf. Part I) 

nr'" +f cl = fr' , (112) 

where / c i (s) is the elastic force acting at arc length s. 
At the ends, where f e \ points in longitudinal direction to 
cancel the external pulling force, one obtains from the 
projection of Eq. (|112J) onto the transverse axis 

r"[ = fr' ± = fr' ± + 0(e) (at the ends) . (113) 

Thus, r'" is typically non-zero at the ends because the 
slope r' — 0{e 1 / 2 ) fluctuates. This, however, corresponds 
to a transverse force at the end, which results in a trans- 
verse deformation. The corresponding bulge of contour is 
only visible on the micro-scale because it spreads with the 
transverse correlation length i±(t) ~ t 1 / 4 . Nevertheless 
this deformation may dominate the growth of the end- 
to-end distance, as is demonstrated for Pulling: With 
the transverse bulk susceptibility scaling as ~ t -1 / 4 
we estimate the magnitude of transverse deformation 
Arj_(s = 0, t) = rj_(0,t) — rj_(0,0) induced by a trans- 
verse force of magnitude ord(f e 1 / 2 ) by 

ord«|Ar x (0, t)|» ~ fe^Hxit) * e 1/2 i 3/4 f . (114) 

The displacement of the end couples to the projected 
length because of the mismatch of the end-tangent 
with the 1 1 -axis by a small angle of typical magnitude 
(l r j_( s = 0)1) — ord(e 1 / 2 ). The expected growth of end- 
to-end distance due to this effect is therefore estimated 
by 

/AJ?J /c J>(t)-(A^||>(t)~ 

(\r' ± (s = 0)|> x (|Arj.(0,t)|> * rt 3/4 f = ^i 3 / 4 f . ^ 

*Cp 

This dominates over the growth law of the bulk, which 
scales like i 7 / 8 / sjlp on short times. The only way to 



avoid the outlined effect is to apply the external force 
strictly tangentially to the end-tangents, which is how- 
ever somewhat unrealistic. The same problem will ex- 
perimentally arise in the short time limit of Release if 
the pre-stretching force was not applied strictly tangen- 
tially. However, the long-time limits arc unaffected by 
this subtlety because the bulk contributions dominate 
over contributions from the end. 

To our knowledge, these end-effects have so far masked 
the subdominant t 7 / 8 -contribution in experiments that 
monitored the time-dependent end-to-end distance (we 
note that in Ref. [2l| the 7/8-scaling is inferred from a 
corresponding scaling of the measured shear modulus of 
an active gel). As we outline in Sec. IVIBI force spec- 
troscopy, on the contrary, may allow to measure the ten- 
sion dynamics, which is itself truly independent of the 
boundary conditions imposed on the transverse displace- 
ments. 



VI. SUGGESTIONS FOR EXPERIMENTS 

While many experiments have been done concerning 
flexible polymers in external force fields [2^, the 
available measurements on driven stiff or pre-stretched 
polymers is not sufficient to verify our predictions. Most 
of these experiments have monitored the transverse and 
longitudinal resp onse on intermediate times where OPT 
is valid (e.g. j24|). Investigations concerning the longitu- 
dinal short-time dynamics are scarce [l(| |2l], |25j . In the 
following, we propose several assays that might be able 
to fill this gap. 

To facilitate the application of our predictions, we rein- 
troduce the parameters k and £ for the following. 



A. Strongly stretched DNA 

The experimental verification of most of our results 
requires stiff polymers with a total length much smaller 
than their persistence length. A remarkable exception is 
Release in the regime of homogeneous tension relaxation, 
which was discussed in Sec. IIV Al For polymers with 
L S> l v , like a typical DNA molecule, this novel regime 
appears if the pre-stretching force obeys 

f»f, p = K ^=^, (116) 

Zp 

where fi p is the Euler buckling force corresponding to 
the buckling length £ p . For forces much larger than the 
cross-over force fi p , the polymer may be considered as 
weakly bending. For those long, but strongly stretched 
polymers our analysis of Release predicts the following. 

After the stretching force has been released the tension 
first propagates through the filament. This regime [TH 
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Il7j ends at the time 

when the tension has to propagated through the filament. 
Here, tr = (^L 2 £ p /(kBT) scales like the longest relax- 
ation time of a Rouse chain with segment length t p . The 

characteristic time t L <C Tr marks the cross-over to a 
regime where the tension profile is roughly parabolic and 
slowly decays according to the power law fit) oc t~ 2 / z . 
This is associated with the projected length i?y growing 
like 

<AJ2,| (t)> = <i?n (t) - fly (0)) = 18 1 / 3 !, ( ±) V3 . 

(118) 

The above analysis strictly holds for the portion of the 
polymer that stays weakly bending. This is not the case 
at the boundaries, for which we refer to existing theories. 
According to the stem-flower model of Refs. [2(| 113, 
the boundaries will develop in time t a "flower" of arc 
length 

( t \ 1/2 

Slower (*)-£(— J , (H9) 

thereby reducing the end-to-end distance by an amount 
of the same order of magnitude as ^Flower itself. It is 
seen, that for t <C tr the shrinkage of the end-to-end 
distance due to the flower is much smaller than that due 
to the weakly bending part of the polymer (the stem), 
Eq. jni- 

Thus, the evolution of the end-to-end distance should 
be described by Eq. (|118f) even for flexible polymers if 
the pre-stretching force is large enough. Since DNA can 
be stretched by very large forces without un-zipping or 
destroying the covalent bonds, we think that the scaling 
(AR(t)) oc i 1 / 3 should be visible in a i?eZease-experiment 
with DNA. The relevant quantities in an experiment with 
A-phage DNA (as in Ref. |2^) in aqueous solution would 
be 

£ p r* 50nm 
L « 20/xm 

a « 2/im (thickness) 
C_l ~ 4?rr//log(L/a) w 1.3 x l(T 3 Pa s 
U P « 0.08 P N 

TR M 7S . 

Though it might require extreme conditions to reach the 
asymptotic limit in Eq. I|118fl . any experiment with finite 
pre-stretching forces larger than f c should be suitable for 
deciding the question whether or not the stem dominates 
the retraction. One would then compare the data with a 
numerical solution |29j| of the governing equation, Eq. l[T]l. 
which describes the stress relaxation in stiff worm-like 
chains. 
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FIG. 7: Two possible realizations of Towing. In (a) the poly- 
mer's left end is being dragged through the solvent with con- 
stant velocity v, whereas in (b) the optical tweezer is im- 
mobile while the solvent flows with constant velocity v. In 
both experiments, the length £u (i) of the boundary layer is 
derived from the pulling force of the tweezer, which can be 
inferred from the displacement of the end within the focus of 
the tweezer. 



B. Single molecule force spectroscopy 

1 . Towing 

As we already mentioned Sec. IIII CI Towing offers the 
possibility to measure the tension propagation by a force 
measurement. The experimental idea is illustrated in 
Fig. E| (a). The polymer's left end is suddenly pulled, e.g., 
by an optical tweezer whose focus moves with constant 
velocity. As a consequence the pulled end follows the 
laser beam with (almost |38|l constant velocity. By mea- 
suring the deflection of the trapped end from the center of 
the beam one can, in principle, extract the pulling force 
f(t) and thus the size £\\(t) ~ f(t)/((v) of the bound- 
ary layer. However, since the laser beam is moving, it 
might represent some problems to dynamically extract 
the deflection. 

A solution to the latter problem is suggested by the 
following Gcdankcn experiment. Consider the above re- 
alization of Towing in the coordinate frame co-moving 
with the left tip, as illustrated in Fig. |7Jb). Then, it 
seems as if the bulk of the polymer was dragged by a 
homogeneous force field to the right while the left end 
is held fixed by the optical tweezer, just as if the optical 
tweezer was spatially fixed while the solvent was homoge- 
neously flowing to the right. In fact, from the polymer's 
perspective there is no difference between both experi- 
ments. Thus, we propose to graft one end of a polymer 
by a tweezer or by the cantilever of an atomic force mi- 
croscope. Then, a homogeneous force field (electric field 
or fluid flow) is suddenly turned on that pulls the bulk of 
the polymer to the right. The deflection of the tip gives 
the pulling force and hence the length of the boundary 
layer. 



TABLE V: Threshold values I c = V K C± f° r diverse biopoly- 
mers. 



Polymer 


Ic 


microtubuli 


1.6 x 10 2 pNms 1/2 


F-actin 


10 pNms 1/2 


Intermediate fil.s 


4 pNms 1/2 


DNA 


0.9 pN ms 1/2 


2. Onset of the nonlinear regime 



One would like to estimate typical time sales for the di- 
verse regimes of tension propagation and relaxation. As 
we have seen above, most of the time-scales crucially de- 
pend on the total length of the polymer, e.g., £\\(t) oc L s 
in the linear regime. This high tunability is of experi- 
mental advantage because one can adjust the setup to 
the observable time scales. On the other hand, it forbids 
to give typical time-scales for those quantities. 

However, there are characteristic quantities that do not 
depend on the length of the polymer. The probably most 
interesting one is the threshold to the nonlinear regime: 
when the product of driving force and the square root of 
the applied time exceeds a certain value I c the polymer 
response becomes nonlinear, 

f Vt > I c = y/^Cl = y/k B Tt p U ■ (121) 

With persistence len gth s of l v = 7mm, 17/xm, 2/zm, 
50nm for microtubu li I30L l3ll| , F-actin [24| , intermediate 
filaments and DNA respectively, and corresponding 
friction coefficients £_l ~ inr]/ ln(£ p /a) (transverse fric- 
tion coefficient per length of a rod of length £ p ), we have 
evaluated I c for some common biopolymers, see Tab. 
Those values can be used to decide whether the response 
of a given biopolymer under a "typical" time-dependent 
external longitudinal force is predominantly nonlinear or 
linear. During a power stroke, for instance, the molecu- 
lar motor myosin exerts a force of about 5pN on an actin 
filament during a time of roughly 1ms (cycle time of the 
power stroke). Hence, the impulse of about 0.5/ c is some- 
what smaller than I c for actin, so that the actin response 
should be linear. 



VII. SUMMARY 

In this paper, we have studied the tension dynamics of 
a weakly bending semiflcxible polymer in a viscous fluid 
theoretically. Starting from the coarse-grained equation 
of motion for the tension, Eq. Q, we elaborated the 
non-linear longitudinal dynamic response to various ex- 
ternal perturbations (mechanical excitations, hydrody- 
namic flows, electrical fields . . . ) that can be represented 
as sudden changes of boundary conditions. For the vari- 
ous scenarios we identified two-parameter scaling forms, 
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that capture the crossover from linear to nonlinear ten- 
sion dynamics. In the limit of large and small arguments, 
where the equilibrium structure of the polymer is self- 
affine, they were shown to reduce to one-parameter scal- 
ing forms, which could be calculated analytically in most 
cases. The growth law £\\(t) ~ t z of the tension profiles 
could be inferred from the scaling variables of the re- 
spective scenarios. This enabled us to develop a unified 
theory of tension propagation. Not only docs it contain 
all cases (correctly) studied in the literature so far. It 
also identifies their ranges of validity and provides new 
predictions. The recovered known results and our new 
predictions are summarized in Figs. H3 and Tabs. HI ITT1 
and GUI 

Various dynamic regimes should be well realizable for 
certain biopolymers. A novel regime of homogeneous ten- 
sion relaxation is a particularly remarkable result from 
the experimental point of view (Sec. IVI Af) . In con- 
trast to previous expectations, this new regime is pre- 
dicted to dominate the relaxation of strongly stretched 
DNA. Moreover, it is an intriguing question, whether the 
tension propagation laws ^ii (t) govern mechanical signal 
transduction through the cytoskeleton [32l l3^ | . We ex- 
pect that the force spectroscopical methods, proposed in 
Sec. IV1BI might be helpful to answer these questions. 

Inclusion of hydrodynamic interactions merely produce 
logarithmic corrections but would give rise to more in- 
teresting effects for polymerized membranes to which 
our discussion could be generalized with otherwise lit- 
tle change. Other natural generalizations including the 
transverse nonlinear response of polymers Q, quenches 
in the persistence length [T3 | and more complex force 
protocols [34| are currently also under investigation. 

VIII. ACKNOWLEDGMENTS 

It is a pleasure to acknowledge helpful conversations 
with Benedikt Obermayer, who, in addition, corrected 
several pre-factors. This research was supported by the 
German Academic Exchange Service (DAAD) through a 
fellowship within the Postdoc-Program (OH) and by the 
Deutsche Forschungsgemeinschaft through grant no. Ha 
5163/1 (OH) and SFB 486 (EF). 

APPENDIX A: CROSSOVER SCALING 
(DETAILS) 

1. Deterministic relaxation on long times (t ^> if ) 

We have for c = 
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and for c = 1 



A- l 'S 1 



3 2^ 9 2 + 1 



dq 1 



D 27T 9 2 + 1 
if d>- 2 



= o(t 



(A2) 



As indicated, both expressions, Eqs. IjAll IA2|) . go to 

zero for large r if 0~ 2 = o(r _1 ). The latter, however, 
follows from the assumptions stated in the main text, 
namely that the tension satisfies the scaling form Eq. I|28|) 
with the requirement in Eq. (|27|l . N.b. it turns out that 
<j> = 0{t) for Pulling and Release and <f> = 0(t 4 ' 3 ) for 
Towing . 



2. Thermal excitation on long times (t ^> tf ) 

We want to show that it is justified to calculate (the 
negative of) the thermally generated stored length rep- 
resented by the term B, Eq. I|47bfl . on long times quasi- 
statically for the semi-infinite polymer. To this end, we 
first insert the scaling ansatz, Eq. for 0(<t,t), 

b=-a r d JL q > fdr 

7a-i2tt J (A3) 

x e -2 9 2 [ 9 2 (r-f)+r° + 1 (0(O-(f/r)° + 1 0K(r/f)--])] 



TABLE VI: Some important notations 


Symbol(s) 


General Meaning 


L 


total length of the worm-like chain 




n 


bending stiffness 




e 


small parameter, defined such that = 


0(6) 


r±{s,t) 


transverse displacement; r = C(e 1//2 ) 




r||(fl,t) 


longitudinal displacement; r\\ = 0(e) 




g(s,t) 


stored- length density; g = r'l/2 + 0(e 2 ) 


= 0(e) 


R 


end-to-end distance 




R \\ 


end-to-end vector projected onto the Ion 


5. axis 




equal up to numerical factors of order 1 




oc 


proportional to 
asymptotically equal 




/ 


persistence length 




k B T 


thermal energy 




f(s,t) 


line tension 




e ± (t) 


equilibration scale for transverse bending modes 


hit) 


scale of tension variations 




f 


external force 




6(s,t) 


thermal force per arc length 





where we introduced the scaling variable £ = cr/r z . Then 
we substitute r — > xr and q — > qr^ 1 / 4 , 



B = -At 1 ' 4 



tlx 



— q 2 

/ A -1 t 1/42TT 
x e -2, 2 [ 9 2 (l-l)+T° + 1 /2(^ ) _ 1 » + l^ ((l - ) )] 



(A4) 



Note that for a > —1/2 (as assumed in Eq. (|27fl 'l the 
factor t q+1 / 2 in the exponent diverges in the long-time 
limit. Hence, for any given wave number the x-integral 
will be dominated by x close to 1 for large enough r ^ 1 . 
This allows us to linearize the exponent in 1 — x when 
performing the x-integral for this given wave number. 
In contrast, for a given time r and 0(£) = 0(1), the 
exponent can be linearized only for large enough wave 
numbers, q ^ g*(r) = r - Q / 2 - 1 / 4 j for which the factor 
q2 T a+i/2 ^ 1 \ n t ne exponent of Eq. i|A4(l is much larger 
than 1. 

Since the integral over q runs over all g-vector, we have 
also to care about the small wave numbers, for which the 
exponent cannot be linearized in 1 — x. To this end, we 
split the g-integral at a wave vector K satisfying 



ql /3 » K » q* 



-a/2-1/4 



(A5) 



which can be found in the limit t>1 under the premise 
of Eq. J23), a > -1/2. The first inequality in Eq. l|A"5fl 

is required for reasons that become clear later on. For 
the upper part B > of the integral we can linearize the 
exponent in 1 — x, 



A' 



2tt 



dx 



x e 



-2q 2 [q 2 (l-x)+T a + 1 ^{^)-x°' + 1 i{i,x-'))] 
-2, 2 (l-x)[, 2 +r 1 / 2 9 T (r"+ 1 ^a/T=))] 



_ 2r i/ 4 f°°dq i^ e -^ 2 b 2 +- 1/2 9r(r a+1 k^))} 
Jk 2tt : 



(tr/r«) 



(A6) 



where we eliminated the scaling variable £ = a /t z , again. 
Using the second inequality in Eq. 1|A5[) it is seen that 
we can drop the exponential for the "interesting" regime 
(f>{a/T z ) = 0(1), where the tension has an appreciable 
value. Inserting back <p = Ta+1 4>{v, T ) we obtain 



B-~ 



-2T 1 / 4 



J dq 



1 



K 2ir q 2 + T 1 / 2 d T <t>(cr,T) 
r°° dq 



1 

~2 x /9 r 0(a,r) ' 



(A7) 
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To obtain the last asymptotics, we have approximated 
the lower bound of the integral by zero. This can be 
justified by the first inequality in Eq. llA5(l . 



K 



= KO( T -*-^ = g*) < qt /3 <C 1 



The remaining lower part B < of the g-integral in 
Eq. (| A4|> is estimated to be small as compared to £?>, 



K dq 



B < = -At 1 '* 

x e -2 q 2 [ q 2 (l-x)+r" + 1 ^(4>(i)-x c ' + 1 ^x^))) 



I ^ q / dX 



27r Jo 



= -^ T l/4^3 /(27r) 



(A8) 



where the first inequality in Eq. (IA5|) has been applied. 
Therefore. B is asymptotically given by Eq. i|A7(l . 



APPENDIX B: RELAXATION OF A 
COMPLETELY STRETCHED POLYMER 

In this section, we consider more closely the interme- 
diate asymptotics 



/0C (i7 



that is approached in the limit 



2/3 



4 7- 8 



7 f 3 / 2 
l p.I< 



C 4 i 



(Bl) 



(B2) 



of Release, which we analyzed in terms of a quasi-static 
approximation in Sec. IV Bl The purpose of this section 
is to justify the quasi-static assumption in the limit of 
/< — > oo (i.e. £f — * 0) where we start with a completely 
stretched polymer and all stored length is generated by 
the action of stochastic forces. 

To this end, we show that, in the limit t — * 0, the 
change in stored length ( A~q) given by Eq. © for the force 
history given by Eq. i|Bl \ asymptotically approaches the 
value one obtains from the quasi-static calculation, 

(^)w=r^{^(^ 292[92t+F(t)l -i) 

+2q 2 /^e-V^Ct-t^+FW-F©] I 

'<^°° 2 q 2 J o * d f e -2 9 2 [ 9 2 (*-*)+F( t )-f (t)] (B3) 



dq 



TT^p 9 2 +/(t) 2l p 



h [/(*)] 



-1/2 



for t -> 



This will comprise an a posteriori justification of the 
quasi-static assumption that entered in the derivation of 
the right-hand-side of Eq. f7H|l . 

The argument closely follows App. IA 21 Inserting the 
force history 



F(t) = / dif(i) = ct a+1 

Jo 



with 



a = -2/3 < -1/2 and C 



Id 



2/3 



(B4) 



(B5) 



into Eq. 1|B3|) and changing variables t — > xt and q — > 
gt -1 / 4 yields 

(Ag)(t) = 2t^ [°° ^-q 2 

x f 1 d ie -^[' > M+^ 1 "(i-^)] (B6) 



As in Sec. IA 21 an approximation to the integral can be 
found in the limit 



Ct a+1 ' 2 » 1 . 
We split the g-integral at K satisfying 



(? y 3 »A'»^(cr+ i / 2 )- 1 / 2 



(B7) 



(B8) 



which can be found in the limit t <C I (=> g* <C 1) 
because a < —1/2. The upper part of the integral 

/•OO 

<AB>)(t) = (...) / dg(...) (B9) 

is dominated by values of x close to 1 and we can linearize 
the exponent in 1 — x, 

(Ag>) (t) = 2t^ 4 r^-q 2 Cdx 
Jk tt Jo 

x e -2g 2 (l-z)[? 2 + (a+l)Cr + 1 / 2 ] 

1/4 f°°dq 1 „ e -^ 2 (l-x)[l 2 + ( a +DCt^] 



= T l/4 



K 



n q 2 + (a + l)Ct a +V 2 

dq 



1 



,1/4 / ^_ 

J K 7T q 2 + (a + \)Ct a + l l 2 

1 r dq 1 



W( a + l ) Cta JK/^/ia+iw+mln T + 1 

1 

~ 2/ p yM 

(BIO) 

where the asymptotics follow from both inequalities in 
Eqs. ifETl lB"8)l. The lower part (A^) of the integral is 
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estimated to be subdominant as compared to (Afp 



(Ag<)(t) = (...) / dq(...) 



K 



< 



2^/4 f K dq 



P Jo 



■q 2 I dx 



3irl p 
oc (Ag>) (t) , 



(Bll) 



crossover at can be understood as follows. For t > t t 
the "speed" of the structural relaxation is determined 
solely by the relaxation times of the bending modes, 
which are related to the transverse friction while the lon- 
gitudinal friction is irrelevant. In contrast, for t <C 
the longitudinal friction substantially limits the speed of 
the relaxation. This suggests to estimate the time 
as follows. From the continuity equation, Eq. (Q, de- 
rived via the multiple-scale perturbation theory (MSPT), 
we can estimate the order of magnitude of the correc- 
tion 5f(s,t) = f(s,t) — / OPT to the flat tension profile, 
Eq. by 



where the first inequality in Eq. I|B8|) has been applied. 
Therefore (A7j) (f) is asymptotically given by Eq. l|B10)l . 

Finally, we want to emphasize the central condition for 
the validity of the quasi-static approximation, 



Cta + l/2 ^ ( & 



2/3 



,-1/6 



» 1 



(B12) 



5f*Cd t (Ag) (/ OPT ,t)L 



(CI) 



OPT can only be applicable if the correction Sf has neg- 
ligible effect on the evolution of the stored length, 



or t < U with U = {L 2 /£ p ) i as given by Eq. 1)73(1. 



(Ag)(f OPT + 5f,t) 



(Ag) (/O pT ,t) 



< 1 , for t > U . (C2) 



APPENDIX C: DEFINING U 

As anticipated Part I, there is a problem-specific time 
t+ limiting the short-time validity of OPT. Physically, the 



The time for which the left hand side of Eq. (|C2)I becomes 
of order may thus be identified with the time t+ before 
which OPT is not valid. 
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